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ABSTRACT 

We present a comprehensive flux resolved spectral analysis of the bright Narrow line 
Seyfert I AGNs, Mrk 335 and Ark 564 using observations by XMM-Newton satellite. 
The mean and the flux resolved spectra are fitted by an empirical model consisting of 
two Comptonization components, one for the low energy soft excess and the other for 
the high energy power-law. A broad Iron line and a couple of low energies edges are 
required to explain the spectra. For Mrk 335, the 0.3 - 10 keV luminosity relative to 
the Eddington value, Lx/L^dd, varied from 0.002 to 0.06. The index variation can be 
empirically described as F = 0.6 logio ^x/^Edd + 3.0 for 0.005 < Lx/Lsdd < 0.04. At 
Lx/LEdd 0.04 the spectral index changes and then continues to follow F = 0.6 logio 
Lx/L Edd + 2.7, i.e. on a parallel track. We confirm that the result is independent of the 
specific spectral model used by fitting the data in the 3-10 keV band by only a power- 
law and an Iron line. For Ark 564, the index variation can be empirically described as 
F = 0.2 logio Lx ji^Edd + 2.7 with a significantly large scatter as compared to Mrk 335. 
Our results indicate that for Mrk 335, there may be accretion disk geometry changes 
which lead to different parallel tracks. These changes could be related to structural 
changes in the corona or enhanced reflection at high flux levels. There does not seem 
to be any homogeneous or universal relationship for the X-ray index and luminosity 
for different AGNs or even for the same AGN. 
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1 INTRODUCTION 

Active Galactic Nuclei (AGN) are known to emit X-rays as 
a result of physical processes active in the innermost regions 
near the central super-massive black hole. X-ray spectra of 
a typical AGN in the 2-10 keV range show primarily the 
signature of a power-law continuum and an iron line. At 
energies < 2 keV there is often a soft excess over the power- 
law emission. The power-law emission is widely considered as 
the outcome of inverse Gompton scattering of thermally pro¬ 
duced accretion disk seed optical /UV photons by a corona 
of hot electrons close to the disc ([Haardt fc Ma raschilllOTU . 

1 19931 : IZdziarski. Poutanen. fc Johnsonll200d ~l. However, the 
geometry and size of the corona as well as the physical mech¬ 
anism governing the energy transfer between the two phases 
are not well understood. 

While, the average value of photon index of the pri¬ 
mary power-law for AGN has been found to be F ~ 
1.9, there is a large variat i on with F rangi ng from 1.5- 
2.5 llNandra fc Pound3ll994l : IPaee et For a given 
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AGN, F varies as a function of time. Understanding the be¬ 
haviour of F with other source properties is expected to give 
important insight into the nature of the radiative mechanism 
and the geometry of the inner regions. 


Previous studies have shown the existence of a corre¬ 
lation bet ween the X-ray pho t on index F and the source 


flux e.g.. llPerola et al. 

Il986l: Singh. Rao & Vahial 

1991; 

Mushotzkv et al.l 1 19931 

iDone. Madeiski. & Zvckil 

200C; 

Nandra & PaoadakisIbOOlll. For many Seyfert 1 AGNs, the 


power-law index shows significant variation and generally 
follows the trend of steeper F with increasing source inten¬ 
sity. During 1997, MCG 6-30-15 was observed with RXTE 
for a duration of 8 days ~ 910 ks. The range of F ob¬ 
served was 1.8-2.2 with flux variation of ~ 3.3 in the en¬ 
ergy band 2-10 keV. This observation showed a tight corre¬ 
lation between the photon index F and the source flux in the 
sense that the power-law st eepens as the source gets brighter 
ll Vaughan fc Edelsonl200lll . The steepening of F (range 1.89- 
2.34) with flux has also been observed in IRAS 13224-3809, 
where a 10 day long ASCA observation revealed a change in 
the power-law flux by a factor of ~ 3.2 in the 2-10 keV range 
llPewangan et al1l2002l l. RXTE observed NGG 5548 on five 
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occasions in 1998 during which a clear positive correlation 
between the photon index (range 1.75-1.93)an d the 2-10 keV 
flux was seen with flux variation of ~ 1.84. (|Chiang et ^ 
l200dl . For the soft energy band (0.3-2.0 keV), similar trend 
between photon index (range 1.7-2.0) and flux was observed 
for a 2 day long BEPPOSAX observation of 3C 120 with flux 
variation of ~ 1.5 (IZdziarski fc Grandil 1200 il l. The Seyfert 
galaxy NGC 4151 was observed on seven occasions by Ginga 
during May 1987-January 1989 (~21 months). In this obser¬ 
vation significant flux variation (~ 3) has been seen along 
with the correlation betw een the photon index and flux 
llYaaoob fc Warwicklll99lh . The 2009 XMM-Newton obser¬ 
vation of Mrk 335 has also shown ‘softer when brighter’ 
property (iGrupe et aLllioi^ l. In all the above examples of 
AGNs, the flux variation during the observations does not 
exceed a factor of ~ 4. 

Studies of a sample of AGNs have also provided ev¬ 
idence that there is significant positive correlation be¬ 
tween the X-ray photon index and the bolometric Ed¬ 
dington ratio, A = LBoi/pEdd , whe r e Lsdd i s the 
Eddington luminosity (iLaor et alJ ll99^;_llm^_J^ lW99|; 
IWang. Watarai. fc Mineshigel 2004 Shemmer et al.l [2006li . 
With larger samples involving sources with higher redshifts 
and greater luminosities, it was observed that the hard X-ray 
photon index cor relates with the bolometric Eddington ratio 
when A > 0.01 teormMLet_aLj|2Q0j; [Shemmeret_aL 2001; 


Saez_et_^ 2008|JS^olewsk^^_Papadnkyj2009||_lC^ 200£ 


Zhou. Xin. fc Zhaol 2010l : Liu et alT 20121 : iBrightman et alJ 

2013h but when A < 0.01, especially for low -luminosity 


AGNs, anti-correlation is seen between T and A (lGu__fc_Cao| 
I2OO9I : [^nstantin et al.l[2009l : lYounes et al.ll201ll : Xul 201 ill 


In general, there have been several studies to investi¬ 
gate the relation between the p hoton index and A for 
differ e nt samples of AGN (e.g . iRisaliti. Young, fc Elvis! 
I 2 OO 9 I: [j in. Ward, fc Don3 12012 1. Studies incorporating 
ROSAT and ASCA observations showed the correlation 
between T and full width at half maximum (FWHM 
of the H/1 emission line dBoller. Brandt, fc Fin i^ll99' 


iBrandt. Mathur. fc Elvisll 19971 : Dewangan et al. 2002l l. 

There is a need to study the relation between spec¬ 
tral index and Eddington ratio spanning a larger range in 
L/pEdd- This is provided by XMM-Newton data for the 
bright and highly variable Narrow Line Seyfert 1 (NLSl) 
galaxy Mrk 335. For comparison we also analyse the exten¬ 
sive data available for another NLSl, Ark 564. 

Narrow line Seyfert 1 galaxies (NLSl) form a sub¬ 
set of AGN which exhibit exceptional features in terms of 
emission-line and continuum properties. Unlike Seyfert 2 
galaxies which show narrow optical emission lines, NLSl 
show the broad emission-line optical spectra of Seyfert 1 
galaxies, but with the narrowest Balmer lines from the 
broad line region (full width at half maximum (FWHM) < 
2000 km s“^ with relatively wea k [OIII]A5007 emission) and 
prominent optical Fe II emission ([Oster brock fc Poggell 19851 : 
IVeron-Gettv. Veron. fc Goncalved 2001^ . NLSl also show a 
number of other extreme properties in X-rays e.g., strong 
soft excess emission below 1 keV, steep 2-10 keV power- 
law continuum and very rapid and large X-ray variability 
(lLeighlvlll999l : [^llo et al.ll2004l ~). Their X -ray spectra often 
present complex behaviour with the presence of cold and 
ionised absorption, pa rtial covering and reflection compo¬ 
nents ( jKomossal 120081 ') and references therein). Also NLSl 


galaxies as a class follow the Mbh —<y» relation if the widths 
of emission lines not strongly affect ed by outflow compo - 
nents are used as a surrogate for ct* (|Komossa fc Xd|2007l 'l. 
Due to all these extreme and ambiguous properties, NLSl as 
a special class of AGN seem to challenge the Unified model 
and need a more careful investigation. Recent observations 
and surveys have pointed towards th e presence of smaller 
black hole masses in N LSl galaxies dBarth. Greene, fc Hoi 
I 2 OO 5 I : iBotte et al.ll20n3 l and which could possibly represent 
an important connection with the less explored intermedi¬ 
ate mass black holes. It is believed that the extreme X- 
ray and other observed properties of NLSl galaxies may 
be due to an extreme value of a fundamental physical pa¬ 
rameter related to the accretion process. This fundamental 
parameter is most likely to be the accretion rate relative 
to the Eddington rate. This is well supported by more re¬ 
cent studies, theoretical considerations and by black hole 
mass estimates from optical emission-line and continuum 
measurements that NLSl are accreting clo s e to their Ed- 


dington rates llBoller. 

Brandtj & FinkI 19961: BorosonI 20021: 

Xu et al.l I 2 OO 3 I: Grune 

I 2 OO 4 I: Iwarner. Hamann. & Dietrichl 

2004lGollin et al.ll2006 

) and hence should be considered im- 


portant testbeds of accretion models. Recent X-ray study of 
NLSl using an optically selected SDSS sample have shown 
that some NLSl show steep X-ray spectra and strong Fe 
II emission while some do not. In this study, a strong cor¬ 
relation was also found between F and the luminosity at 
1 keV, Likev suggesting differences in Ltoi/LEdd among 
the NLSls in the sample (I Willi arns. Pogge. fc Mathuil[2002 : 

IWilliams. Mathur. fc Poggell2004h . 

Mrk 335, also known as PG003-I-199 is a n earby N LSl 
galaxy at a redshift z=0.026 dLonginotti et aHl2007al l and 
has a well measured blac k hole mass of 1.4 x l O^ M(?) from 
rever beration mapping (|Peterson et al.l l2004l : iGrier et al.l 
I2OI2II . It has been the target of most X-ray observatories. 
XMM-Newton observed Mrk 335 for the first time in 2000. 
Analysing the RGS data, an absorption edge at 0.54 keV due 
to Galactic oxygen was reported and the soft excess was de¬ 
scribed as a combination of bremsstrahlun g emission and 
ionise d reflection from the accretion disk (|Gondoin et al.l 
I2OO2II . The XMM-Newton observation was later reanalysed 
and a narrow absorptio n feature at 5.9 keV was found 
l|Longinotti et all l2007al L In 2006 January, XMM-Newton 
re-observed Mrk 335 for 133 ks which revealed a double- 

J leaked Fe emission f eature with peaks a t 6.4 and 7.0 keV 
O’Neill et al.ll2007h . (|Larsson et al.l[2008l ') studied Mrk 335 
using a 151 ks Suzaku observation performed in 2006. They 
modelled the data using a power law and two reflectors in 
which an ionised, heavily blurred, inner reflector produces 
most of the soft excess, while an almost neutral outer reflec¬ 
tor (outside ~ 40rg) produces most of the Fe line emission. 
They also verified their model using the 2006 XMM-Newton 
data and did not see any correlation between photon in¬ 
dex and power law flux. But subsequently a marginal trend, 
i.e., the source becomes softer with in creasing count rat e 
has been reported for the same data dGrupe et aLllioi^ i. 
When observed by XMM-Newton in July 2007 for 22 ks, 
Mrk 335 was found to be in extremely low X-ray flux state 
dGrupe et al.l [200^ . The spectrum of this low flux state 
of Mrk 335 was explained by partial covering and blurred 
reflection models. Mrk 335 was again observed by XMM- 
Newton in 2009 for 200 ks spread over two consecutive or- 
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bits. The X-ray continuum and timing propertie s have been 
descr ibed by using a blurred reflection model (iGallo et al.l 
12011^1. The 2009 XMM-Newton data have also been analysed 
bv lcrupe et al.l (l2012ll along with the Swift data. Partial ab¬ 
sorption and blurred reflection models gives equally good fit 
to the spectrum. The number of observations of Mrk 335 by 
XMM-Newton over time provides an opportunity to study 
the variation of the high energy photon index with luminos¬ 
ity in a systematic manner. _ 

At a redshift z=0.02469 l|Huchra et al.l[l999l '>. Ark 564 
is the brightest NLSl galaxy in the 2.0-10.0 keV 
range , L(^ 2 -io)keV = 2.4 x 10"^® erg s“^ dTurner et al.l 
I 2 OOII). Ark 564 has been studied acros s all wavebands 
llShemmer et al ] |200ll : [Romano et al.ll2003l . In the 2000 and 
2001 XMM-Newton observations of Ark 564, IVignali et al.l 
l|2004h reported an edge-like feature in the EPIC data at 
~ 0.73 keV and interpreted it as the OVII K absorp¬ 
tion edge. The ~100 ks 2005 XMM-Newton observation of 
Ark 564 have been described as either a power law and two 
black-bodies or a relativistically blurred p hoto ionised disk 
reflect ion model (jPanadakis et al] l2007l l. IPewanean et al.) 
ll2007^ analysed the 2005 data and found two warm absorber 
phases. During 2011, Ark 5 64 has been observ ed by XMM- 
Newton for eight occasions. iLegg et al.l (|2012l l confirmed a 
significant soft lag in the 0.3-1.0 keV and 4.0-7.5 keV bands 
and suggested a distant reflection origin. Thus, Ark 564 has 
been also observed several times by XMM-Newton making 
it another good candidate. 

Here we investigate the dependence of power-law index 
with the X-ray Eddington ratio for the combined EPIC spec¬ 
tra of Mrk 335 and Ark 564. Using flux resolved spectroscopy 
we study the power law index variation against X-ray Ed¬ 
dington ratio. The paper is organised as follows. Section 2 
describes the observations and the data reduction procedure 
and Section 3 explains the spectral model used for the anal¬ 
ysis. The details of the flux resolved spectroscopy are given 
in section 4. Finally, we conclude this work with a summary 
of the results and discussion in section 5. 


2 OBSERVATIONS AND DATA REDUCTION 

For the analysis, we use all the available archival data of 
Mrk 335 and Ark 564 from the XMM-Newton observa¬ 
tory ([.lansen et al.ll200lh . The list of observations are shown 
in Table 1. The XMM-Newton data have been processed 
in the standard way using the SAS version 12.0. For all 
cases we have consi dered data from the EPIC-pn camera 
llStruder et aLlliooIh only. 

The data were cleaned for high background flares and 
were selected using the conditions PATTERN ^ 4 and 
FLAG == 0. We have checked for pile-up in all cases using 
SAS task epatplot and found that for both the sources some 
observations are affected by pile-up (Table 1). We reduced 
the pile-up by excluding the innermos t source emission using 
an annular region as considered by (iLegg et al.ll2012f l. Ex¬ 
cept for the pile-up affected observations, the source spectra 
have been extracted from circular regions of radius 35 arc 
s centred on the maximum source emission. The redistribu¬ 
tion matrices and auxiliary response files were created by 
the SAS task especget. Spectra were grouped such that each 
bin contained at least 30 counts. Spectral fits to the data 


Table 1. Observation log of Mrk 335 and Ark564 by XMM- 
Newton 


Mrk 335 

OBS. ID 

Duration(s) 

Date 

Pile-up 

0101040101 

36910 

2000-12-25 

yes 

0306870101 

133251 

2006-01-03 

no 

0510010701 

22580 

2007-07-10 

no 

0600540601 

132315 

2009-06-11 

no 

0600540501 

82615 

2009-06-13 

no 

Ark 564 

0006810101 

34466 

2000-06-17 

no 

0006810301 

16211 

2001-06-09 

no 

0206400101 

101774 

2005-01-05 

no 

0670130201 

59500 

2011-05-24 

yes 

0670130301 

55900 

2011-05-30 

no 

0670130401 

63620 

2011-06-05 

no 

0670130501 

67300 

2011-06-11 

yes 

0670130601 

60900 

2011-06-17 

no 

0670130701 

64420 

2011-06-25 

no 

0670130801 

58200 

2011-06-29 

yes 

0670130901 

55900 

2011-07-0 

yes 


were performed with XSPEC version 12.8.0 llArnau dl ll996lL 
The background subtracted light curves have been produced 
with the tool EPICLCCORR which are binned with 400 s 
bins. 


3 PHENOMENOLOGICAL MODEL 

We fit the EPIC-pn spectral data in the 0.3-10.0 keV range, 
by a simple phenomenological model consisting of two ther¬ 
mal Comptonization models. In particular, the soft excess is 
described by the XSPEC model “nthComp” and for the hard 
X-ray e mission, we have us ed XSPEC convolution model 
“Simpl” llSteiner et aLll^OSl l. Galactic absorption is taken 
to be Nir=3.99x 10^° cm“^ for Mrk 335 and Nir=5.34x 10^° 
cm“^ for Ark 564 (iKalberla et al.ir2005l f. For Mrk 335 there 
are clear residuals at around 6 keV signifying the presence 
of a broad Iron line which we modelled using the “diskline” 
llFabian et al.|[l989l L The inner radius of the line emitting 
region Ri„ is fixed at 6rg and the outer radius Tout is fixed 
at lOOOrg. The inclination of the disk is fixed at 40°. For the 
“nthcomp” model we find that the results are not sensitive 
to the input seed photons which we fix to be a blackbody at 
0.05 keV. For both the sources, addition of 2 or 3 low energy 
absorption edges improves the spectral fitting. For one of the 
Mrk 335 observations there is also a hint of a narrow emis¬ 
sion line at 7 keV. Figure [T] presents the unfolded spectra for 
the observations 0101040101, 0306870101, 0600540501 and 
0600540601 of Mrk 335. The best fit spectral parameters and 
the reduced for all the observations are listed in Tables 
[2]and|3l 

The reduced x^ of the fits to the observations range 
from 1 to 1.6, which suggests that the underlying spectra 
may be more complex than the phenomenological model 
used here, especially at low energies. However for the moti¬ 
vation of this work, the phenomenological model used here 
is adequate since even for a reduced x^ ~ l-6i tho model 
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Figure 1. EPIC PN unfolded spectra and residuals (Ax) of Mrk 335 observations i.e. 0101040101 (top right), 0306870101 (top left), 
0600540501 (bottom right) and 0600540601 (bottom left) for the phenomenological model with best fit parameters listed in Table [2] 



Time 





Time Time 


Figure 2. EPIC PN light curve (in 400 s bins) of Mrk 335. The light curves were extracted in the energy range of (0.3-10.0) keV. The 
time ranges for flux resolved spectroscopy have been selected corresponding to different counts s as shown by the horizontal dotted 
lines. 
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Figure 3. EPIC PN light curve (in 400 s bins) of Ark 564. The light curves were extracted in the energy range of (0.3-10.0) keV. The 
time ranges for flux resolved spectroscopy have been selected corresponding to different counts s as shown by the horizontal dotted 
lines. 


Table 2. Spectral parameters for Mrk 335 derived from different XMM-Newton observations in the (0.3-10.0)keV range. 


Model 

parameter 

0101040101 

0306870101 

0510010701 

0600540601 

0600540501 

zedgel 

Ec (keV) 

0.63(froze) 

0 351+0 003 

U.OOi_Q QQ4 

i.U4^_0 035 

u. 'ou_0 003 

U. 'OU_0 004 


T 

u.uy ^_o 031 

n 99 Q + 0-027 
U.zzo_o 000 

0 359+0 091 
u.ouy_0 003 

0 775+0-027 

‘ ' '3-0.042 

n c^c3+0.029 
U.OOO_0 017 

zedge 2 

Ec (keV) 

1 1 97+0.042 

-^•J-^<_0.046 

±.1^U_Q 230 

j..ouy_o 003 

u.y4:o_0 012 

u.y^z_o oil 


T 

U.UOO_Q Q31 

U.UZD_Q QQg 

n 077+0.101 
' -0.095 

n “^1 Q+0.030 
u.o±o_0 041 

n 90Q+0.024 

u.zoo_0 040 

zedge3 

Ec (keV) 

- 

- 

0.686(Froze) 

- 

- 


T 

- 

- 

u.oz^_o 100 

- 

- 

Diskline^ 

E, (keV) 

c o-| 

u.o±o_Q 1^4 

a 97c-t-U.Ucs:fi 
O.Z/0_Q Q41 

o.uuy_0 0gg 

u.u^z_o 061 

o.yoz_o 064 


/3 

— 2 609"^^’^'^^ 

z.Duy_5 gQQ 

1 09^+0.643 

1.9ZD_q 414 

y.bDo_g 00g 

i .^1D_4 763 

— 5 061 

O.Ubi_4 476 


Norm(x 10“^) 

0 44Q+0-254 

U.z± i _o 033 

u.o4y_0 076 

u.z/u_0 033 

u.zyz_0 022 

zgaussian 

Ec 

- 

i •UU^_Q Q31 

- 

- 

- 


Norm(x 10“®) 

- 

U.UU^_Q QQ2 

- 

- 

- 

Simpl 

r 

^•^^^-0.049 

1 oqc+O.OOQ 

i.5»b_Q QQ5 

l.UO^_0 046 

l.D/D_0 007 

1 7 Q 7 ”t“U.UlS 
-0.008 


FracScat 

0 iQfi+O-Oll 

U.10U_Q QQg 

o.ostIHm 

U.O1U_0 104 

n 900+0.012 
u.zoo_0 007 

0 193+0-007 

u.iyo_0 006 

nthComp 

r 

9 0^9 + 0.052 
4.5b4_Q Qgy 

0 q'7't+0.008 

'-0.023 

1 QQn+0-094 

i.y5U_0 002 

1 099+0.087 
i.yzz_o 213 

^.u^o_0 219 


kT(keV) 

U.ZOO_Q Q31 

n Q9Q+o.oi9 
u.o^o_Q Q20 

U.Z1U_0 010 

n 1 cro+0.008 
U.iOO_0 0ig 

0 159+0-008 
u.ioy_o 020 


Norm(xl0“^) 

n 49Q+0.029 

U.42»_q Q32 

U. 037 

u.uyo_0 030 

U.UOU_0 00g 

u.uyi_0 009 

reduced y^/d.o.f. 


1,04/144 

1.48/163 

1.70/130 

1.40/161 

1.52/159 


'f’The values of xin and Tout have been fixed at Qvg and lOOOr^ respectively while the inclination has been fixed at 40°. 
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Table 3. Spectral parameters for Ark 564 derived from different XMM-Newton observations in the (0.3-10.0)keV range. 


Model 

Parameter 

0006810101 

0670130301 

0670130801 

0670130901 

0670130201 

0670130601 

zedgel 

(keV) 

U.OUQ_o.oi5 

n 708+0 008 

U. <^Uo_o.o08 

0 84 Cv +0 013 

U.CS40_o.oi5 

n 790+0.010 

U. <^0_Q QIQ 

0 990+0-011 

U.OOU_Q QJ^g 

0 947+0.005 

U.O4/_0 QQg 


T 

n n 1 +^-^22 
J--Lr_o.022 

n io'^+o oi2 

U.iUO_0.0l0 

0 244+0-050 
U.^44i_o 051 

u.uyo_Q Q13 

u.ooy_o.o44 

u.zDy_o 021 

zedge2 

Ec (keV) 

u. i rz_Q QJ3 

l.UO4_0.045 

0 71 O+0.012 

U. i ro_o.o06 

1 071 +0.053 

i.ui r_o.062 

0 7“^! +0-020 
U./OL_q Q23 

0.714«;““ 


T 

Q -14-7+0 030 
-0.022 

n r)QQ+0 013 

u.uoo_o.oi4 

u.u/z_o.oio 

0.029t°;°i® 

U.U40_Q QQ7 

0.063l°+i2 

zedgeS 

Ec (keV) 

1 1 99+0-126 
L.iZZ_o 094 

-1 QQ-|+0.038 
i.Z9i_o.o39 

l-l60l“;°« 

1 991 +0.159 

L.ZOi_Q Qg7 

1 91 7+0.033 
i-^i<_0.031 

1 1 97+0.021 
l--L^<_0.019 


T 

U.UZZ_Q 022 

U.U4(_o 014 

0.086l°;°J? 

U.UZO_Q Qlg 

0 0^7+0.014 

U.UO/_Q q;L4 


Simpl 

r 


9 c-j o+o.oi4 
z.oiO_o.oi3 

9 41 q+0.014 

4.41»_o 008 

2 40^1+0-018 
z.4yo_Q Q17 

9 zl48+0-0i9 
4.445_Q QQQ 

9 /1 98+0-015 
ZAZQ_o_oi4 


FracScat 

n iqq+0.008 

u.iyo_o.oo7 

0 1Q2+0-004 

U.iy-^_0.oo5 

0 1 e;9+0 004 

U.iOO_o.o04 

0 1Q4+OOO6 

u.iy4_o 006 

n 1 ^9+0.005 

U.10Z_Q QQ5 

^ -,.-,+0.004 

U.iOi_o 004 

nthComp 

r 

9 9QC+0.052 

z.zyo_o.o47 

9 9Q9+0.010 

2 c7«+0-059 

Z.D/D_o.o50 

9 oc:o + 0.016 
Z.OOO_Q Q22 

9 7qc;+0.063 

9 q/:?9+0.051 

4.5bZ_o 044 


kT(keV) 

n 1Q4+0-020 

^•-*-^4-0.020 

0 lQCv+0 009 

u.iyo_o.oo8 

u.zzu_o.o20 

0 1Q2 + 0-011 

u.iyz_o.oio 

U.Z4U_q Q23 

^•^'^^-0.022 


Norm 

0 014+0-0008 
U.Ui4_Q Q009 

0 01 9+0.0003 

u.Uiz_o 0003 

0 011 +0-0004 

u.uii_o.oo03 

0 01 9+0-0004 

U.U10_Q QQQ3 

0 01 [^+0-0005 

U.U10_Q QQQ4 

0 011 +0-0004 
U.Uii_o 0003 

red /d.o.i. 


1.03/149 

1.58/161 

1.62/160 

1.17/158 

1.03/159 

1.30/164 

Model 

Parameter 

0670130501 

0006810301 

0206400101 

0670130701 

0670130401 


zedgel 

Ec (keV) 

0 ‘:i4fi+0 006 

U.04D_Q QQY 

0 Qoo+0.006 

u.ooo_0.007 

n K-| 0+O.OO6 

U.010_Q QQg 

0 [^Q7+0-045 
u.ooi_o 037 

0 S41 +0-008 

U.04i_Q Q^l^g 



r 

U.4UD_o.o44 

0.595t°of2t 

U.UyO_Q QQQ 

u.uoo_Q Q20 

0 909+0-041 

U.OUO_Q Q4J^ 


zedge2 

Ec (keV) 

0.718l°;°J= 

0.489l“;°i2 

0 700+0-005 

U. QQ7 

0 710+0-013 

U. mD_o oil 

0 710+0-009 

U+iU_Q QQ5 



T 

0.069l“;°“ 

0.188lg;°24 

oil 0+0-011 
U.ilD_o.012 

U.iZD_Q Q24 

U.UlO_Q q;L2 


zedgeS 

Ec (keV) 

l-197taoi 

0 710+0 015 

U. mu_o.oi7 

-,-,.+0.018 

i-J--Li_0.019 

1 149+0.023 

l.iho_q Q27 

1.164l0;°|° 



T 

0.068l“;°i2 

U.iU^_0.oi6 

U.U04_Q QQg 

0 Ofil +0-014 
U.UDi_o 014 

0.046i°;°ll 


Simpl 

r 

9 407+0.oil 
Z.^O(_o.oi4 

9 07^+0.069 

z.o/u_o.oi2 

9 479+0.009 

z.H/z_o 004 

2 49C+O OI6 

Z.40D_Q QJg 

2 410+0-012 
^.4iO_Q QJ^5 



FracScat 

n 144+0-004 
U-l^^_0.004 

0.085lg;°i5 

n 1 qo+0.003 

U.iOO_Q QQ3 

0 1 Q9+0.004 

U.iOCS_Q QQ5 

0 1 [^0+0-005 

U.10y_Q QQ4 


nthComp 

r 

9 QQ9+0.029 

z.oyz_Q 029 

Q iq>^+0.024 

9.iZl_o 021 

9 981 +0.019 

z.ooi_o.o22 

9 1 7^+0.022 
z.l/o_q Q27 

9 740+0.049 

4./4«_q Q5j^ 



kT(keV) 

u.zou_o.o23 

U..Z4U_o.oo4 

n 1 88+0-009 

U.iO»_Q QQg 

0 101 +0-009 

U.iyi_o 008 

n 900+O.032 

U.zoo_Q Q22 



Norm 

0 m 2+0 0004 

u.uiz_o 0002 

u.uuy_o.oooi 

0 010"*"0 0001 
U.UiU_Q QQQ3 

U.UUy_Q QQQ3 

0 011 +0-0003 

u.uii_Q 0003 


red /d.o.l. 


1.35/163 

1.16/148 

1.62/165 

1.51/163 

1.49/164 



described the data well at a few percentage level and most 
of the discrepancies are at low energies. The high energy 
photon index is not too sensitive to the actual model used 
as we show later when we fit only the high energy part (3-10 
keV) of the spectra and obtain qualitatively similar results. 
Nevertheless, we caution against over-interpretation of the 
best ht values obtained in these hts, especially for the com¬ 
ponents affecting the low energy part of the spectra. 


4 THE EFFECT OF COMPLEX ABSORPTION 

AND REFLECTION COMPONENT 

The spectra of Mrk 335, may be affected by complex ab¬ 
sorption and there may be a relativistically blurred reflection 
component which may dominate at low and high energies. 
Given the quality of the data many of these complex mod¬ 
els may be degenerate, however in the present context it is 
useful to quantify the effect of these models on the high en¬ 
ergy spectral index. The motivation here is not to obtain 
a physically self consistent model but rather to understand 
their effect of the high energy spectral index and hence in 
the primary result of this work. 

To take into account the possibility of complex ab¬ 
sorption, we include in the phenomenological model, partial 
ionised absorber re presented by the XSPEC routine “zxipcf” 
llMiller et al.ll200^ . The absorption is characterised by three 
parameters namely the column density, the covering fraction 
and the ionisation parameter of the absorber. For the ad¬ 
dition of three parameters the changes in for the four 
high flux observations with IDs 0101040101,0306870101, 


0600540601 and 0600540501 were 12.7, 22.0, 3.1, 9.0 and 
12.0 respectively. More importantly, the best fit high en¬ 
ergy spectral indexes obtained were 2.05 ± 0.05, 1.78 ± 0.06, 
1.75 ± 0.05 and 1.81 ± 0.02. When compared with the high 
energy spectral indexes obtained from the phenomenological 
model, one can see that the errors are larger and any change 
in the index is less than 0.1. In contrast, for the deep low 
flux state, ID 0510010701, the Ax^ = 45 and the spectral 
index obtained was 1.55±0.15 which is significantly different 
from the phenomenological model. 

Next we consider the possible effect of a complex rel¬ 
ativistically blurred reflection component in the spectra. 
Since, the ionised reflection component produces a complex 
soft excess which needs to be modelled with absorption, we 
limit our analysis to energies > 2 keV. This is adequate 
since our interest here is to study the effect of the compo¬ 
nent on the high energy index. Since for the phenomenolog¬ 
ical model, the spectra shows the presence of a narrow and 
broad Iron lines, we consider a power-law and two reflec- 
tion components repr esented by the table model “reflionx” 
(IRoss fc FabianI [ 2 OO 5 II . For one of the the reflection com¬ 
ponent we convolve it using the relativistic blurring model 
“kdblur” where we fix the inner and outer radii, but al¬ 
low for the emissivity index to be free. The best fit power- 
law index for the four high flux state observations with IDs 
0101040101,0306870101, 0600540601 and 0600540501 were 
2.09±0.05, 2.16±0.02, 1.77±0.03 and 1.64±0.05. We note 
that for the observation with ID. 0306870101 the change in 
spectral index as compared to the phenomenological model 
is 0.28 while for the others it is < 0.15. For the deep low 
flux state, ID 0510010701, the spectral index obtained was 
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Figure 4. The high energy photon index, F versus the X-ray Ed¬ 
dington ratio for Mrk 335. The X-ray Eddington ratio is 
where Lx is the unabsorbed luminosity in the 0.3-10 keV range. 
The solid line is a fit to three of the data sets and has a slope 
m = 0.64 it 0.04 and intercept c = 3.08 it 0.08 with a reduced 
= 0.57. The bottom dotted straight line is a fit only to the 
data set (ID0306870101) and has a slope m = 0.65 it 0.04 and 
intercept c = 2.76 it 0.07 with a reduced x^ = 0.76. 


OBS-ID0306870101 i 
OBS-ID0600540601 ; 
OBS-ID0600540501 i 
OBS-ID0101040101 ■ 


Figure 5. Same as Figure |4] except that only the energy range 
3-10 keV is taken into account. The spectra are fitted with a 
power-law and Iron line. The luminosity L[fx corresponds to the 
energy range 3-10 keV. Note the similarity with Figure [4] which 
implies that the qualitative results are not sensitive to the spectral 
model adopted. 


3.00±.07 which is radically different than the value obtained 
earlier. 

Thus, the effect of complex absorption or reflection is 
dramatic for t he deep low state. This is expected based on 
the analysis of iGrupe et al.l ll2008tl . For the other higher flux 
states, complex absorption may change the index by < 0.1 
from those of the phenomenological model. The effect of 
Complex reflection on the index is also modest < 0.15 but 
we note that it may be large ~ 0.28 for the observation with 


ID. 0306870101, which incidentally is also the highest flux 
state in the sample. 

5 FLUX-RESOLVED SPECTROSCOPY 

The light curves for each of the observations of Mrk 335 
and Ark 564 are shown in Figure [5] and Figure 0 which re¬ 
veals that for nearly all observations there is significant flux 
variation in timescales of ~ lO"^ s. To investigate the varia¬ 
tion of the photon index, each observation was split into 3 
or 4 flux levels (marked by horizontal dotted lines in Fig¬ 
ures [2] and 0 and the corresponding spectrum was gener¬ 
ated. We do not consider the 2007 observation of Mrk 335 
(ID0510010701) for the flux resolved spectroscopy since as 
shown in the previous section the spectra is clearly complex 
and the phenomenological model used here is not adequate. 

Each of the flux level spectra was fitted with the same 
phenomenological model used for the average spectra. How¬ 
ever, due to the lower statistics of the flux resolved spec¬ 
tra, there were several parameters which were either not 
constrained or were consistent within error bars to be not 
varying during the observations. Thus, several spectral pa¬ 
rameters were fixed to their best fit values obtained from 
the fitting of the average spectra. These were, the energies 
of the different edges and emission lines, the temperature 
and normalisation of the soft photon Comptonization (i.e. 
kThb and normalisation of the “nthcomp” component) and 
the normalisation and emissivity index of the Iron line (i.e. 
P and normalisation of the “diskline” component). 

For each flux resolved spectra, the unabsorbed flux in 
the 0.3-10 keV band was computed using the XSPEC model 
“cffux”. The fluxes were converted into luminosities, L by 
using the luminosity distances of 103 Mpc for Mrk 335 and 
98.5 Mpc for Ark 564, which are the average values quoted 
in the NED website. For Mrk 335, the black hole mass was 
taken to be 1.4 xlO^ Mq llPeterson et al.ir2004l : [&ier et al.l 
I 2 OI 2 II while f or Ark 564 we adopt a value of 2.61 x 10® Mq 
( Botte et al.ll2003 l to obtain their respective Eddington lu¬ 
minosities LEdd- The best fit parameters for the flux resolved 
spectra are listed in Table [4] 

In this work we consider the X-ray Eddington ratio i.e. 
Lx jLsdd instead of the Bolometric one and discuss the im¬ 
plication of this in the last section. In Figure ID the high 
energy photon index is plotted against the X-ray Edding¬ 
ton ratio for the different flux resolved spectra. Three of 
the observations cover an order of magnitude range in X- 
ray Eddington ratio from 0.004 to 0.04 and show a tight 
correlation between the photon index and X-ray Eddington 
ratio. A straight line fit to these three observation sets gives 
a slope m — 0.64 ± 0.04, intercept c = 3.08 ± 0.08 and re¬ 
duced — 0.57. The observation with the highest X-ray 
Eddington ratio > 0.04, does not follow the linear relation. 
Instead fitting the flux resolved spectra for that observation 
only, one gets a straight line with slope m = 0.65 ± 0.04, 
intercept c = 2.76 ± 0.07 and reduced x^ — 0.76. In other 
words, the data follows a parallel track. Note that the lowest 
luminosity observation {Lx/Lsdd ~ 0.002) was not included 
in the straight line fit. Nevertheless, the fit to the high flux 
data, passes through the lowest flux one. 

To verify that the results do not depend sensitively on 
the phenomenological spectral model used, we repeat the 
flux resolved spectroscopy analysis using only the energy 
band 3-10 keV. We fit this energy band with a power-law and 
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Figure 6. The high energy photon index, P versus the X-ray Ed¬ 
dington ratio for Ark 564. The X-ray Eddington ratio is Lx/^Edd 
where Lx is the unabsorbed luminosity in the 0.3-10 keV range. 
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Figure 7. Summary figure for the high energy photon in¬ 
dex, r versus the X-ray Eddington ratio for both Mrk 335 and 
Ark 564. The X-ray Eddington ratio is Lx/L^dd where Lx is 
the unabsorbed luminosity in the 0.3-10 keV range. For Mrk 335 
two parallel tracks are observed with a possible switching at 
Lx/Lsdd 0.04, while for Ark 564 the correlation is flatter 
and with significantly more scatter. 
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a broad Iron line. Figure [5] shows the results of this analysis 
where the photon index is plotted against Lhx/LE dd, where 
Lhx is computed using the unabsorbed flux in the 3-10 keV 
band. 

Figure [6] shows the results of the flux resolved analy¬ 
sis (for the complete 0.3-10 keV band) for Ark 564. Here, 
despite the larger number of observations, the flux varia¬ 
tion is modest with the X-ray Eddington ratio ranging from 
0.02-0.07. Although there is more scatter than the case of 
Mrk 335, there is a correlation between the spectral index 
and X-ray Eddington ratio. A straight line fit gives a slope 
m = 0.22 ±0.08, intercept c = 2.75 ±0.1 and a large reduced 
= 2.75. From the Figure, it seems that for the flux re¬ 
solved spectra of each individual observation the correlation 
may be better and perhaps the different observations are 
parallelly shifted with respect to each other. However, the 
statistics is not good enough to make any concrete state¬ 
ments. 


6 DISCUSSION 

Figure [3 summarises the results of the work by showing the 
variation of the high energy photon index versus the X-ray 
Eddington ratio for both Mrk 335 and Ark 564. For Mrk 335, 
there are two parallel tracks with the possibility that the 
sources shifts to the lower track for LxjLsdd > 0.04. For 
Ark 564 the points are clustered at a larger F and while 
there is a positive correlation, the slope is flatter and points 
show much more scatter than the points for Mrk 335. 

It should be noted that the luminosity used in this work 
is in the energy range 0.3-10 keV and is not the bolomet- 
ric one and hence these results cannot be compared directly 
with those where the bolometric luminosity or X-ray lumi¬ 
nosity in some other energy range has been used. The con¬ 
version to Bolometric luminosity often involves several un¬ 
certain factors and hence is avoided in this work. Since for 
Mrk 335, we find the the index correlates well with the X-ray 
luminosity for nearly an order of magnitude change in lumi¬ 
nosity, this may indicate that over this range the bolometric 
correction is nearly constant. On the other hand, it may well 
be that for this source the X-ray index correlates better with 
the X-ray luminosity and not with the bolometric one. The 
uncertainties and model dependency of estimating the bolo¬ 
metric luminosity (including the general non-availability of 
simultaneous multi-wavelength data) does not allow for con¬ 
crete statements. For Mrk 335, there is only one observation 
for which the X-ray Eddington ratio is larger than 0.04 and 
the source follows a lower parallel track. It is necessary to 
confirm this behaviour with future observations when the 
source is equally luminous. It might be that there is a real 
transition at the X-ray Eddington ratio ~ 0.04 or that the 
source can actually follow any of the two tracks at the same 
luminosity. 

In the Comptonization context, the high energy power- 
law index is inversely proportional to the Compton Ampli¬ 
fication factor A, which is the ratio of the luminosity of 
Comptonizing cloud, Lc to the input soft photon luminosity 
Linp i.e. A = LcjLinp. The input luminosity Li„p depends 
on the luminosity of the soft photon source and the fraction 
of photons which enter the Comptonizing region. The lat¬ 
ter depends on the accretion geometry of the system. Thus, 
the correlation between the photon index and the observed 
luminosity could occur if A decreases with luminosity or in 
other words Lc varies less rapidly with the observed lumi¬ 
nosity than Linp- The shift to the lower parallel track can 
be explained if there is a decrease in Linp at the X-ray Ed¬ 
dington ratio of ~ 0.04, perhaps caused by a change in the 
fraction of photons entering the Comptonizing region. This 
would mean that the accretion geometry for the two parallel 
tracks are different. 

The high energy photon index could also be affected by 
the presence of a strong reflection component, especially if 
the reflection is from partially ionised matter and/or is rel- 
ativistically blurred such that it has a significant contribu¬ 
tion to the spectra below 10 keV. The reflection component 
would tend to flatten the high energy spectrum and hence 
the correlation may be caused by the reflection component 
decreasing as the source becomes more luminous. Indeed, 
the anomalous low luminosity observation of Mrk 335 of July 
2007, can be modell ed as being domin ated by a blurred re¬ 
flection component llGrupe et al.ll200^ . It may also be that 
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the primary correlation between the index and luminosity 
is due to variation of the Compton Amplification factor, 
while the shift to a different parallel track is due to appear¬ 
ance of a strong reflection component. The appearance of a 
stronger reflection component would also mean that there 
was a change in the accretion disk geometry. Thus, a qual¬ 
itative change in the accretion disk geometry seems to be 
required to explain the parallel tracks seen for Mrk 335. 
Changes in the accretion disk geometry at the same X-ray 
luminosity seems to indicate that the accretion flow is not 
determined only by the local accretion rate but rather it may 
also depend on the previous history of the accretion rate 
variation. Our brief analysis of including a blurred reflec¬ 
tion component suggests that this indeed may be the case. 
It is particularly interesting to note that for the highest flux 
observation (ID 101040101) a reflection component may in¬ 
crease the index by 0.28 while a more modest variation is 
expected for the others (Section 4). If indeed this highest 
flux observation has a stronger reflection component and its 
spectral index is more close to ~ 2 then it would seem that 
the index has a positive correlation with flux but saturates 
to ~ 2 at high Eddington ratios. Interestingly a similar be¬ 
havior is observed in Fig. 5, where the results for a simple 
power law plus Iron line fit to the 3-10 keV energy range 
are shown. However, this can be confirmed only with broad 
band data covering energies > 10 keV where the reflection 
component will be better constrained. 

The results for Ark 564 are significantly different than 
that for Mrk 335, with a flatter correlation and more scat¬ 
ter. This suggests that the behaviour of AGNs is not homo¬ 
geneous and perhaps cannot be understood by studying a 
large sample of them. Long term sensitive and broad band 
monitoring of individual AGNs may be expected to provide 
better insights. 
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Table 4. Spectral parameters for each flux state (for Mrk 335 Sz Ark 564) derived from different XMM-Newton observations in the 
(0.3-10.0)keV range. The full version is available to download online 


Obs ID 

flux state 

zedgel 

zedge2 

zedge3 

Simpl 

nthComp 

red x^/d.o.f. 



T 

r 

r 

r 

FracScat 

r 


0101040101 

1 

0.002 X 10-4 

U.UM_0.o65 

- 

2 1-14+0.157 

0 14Q+0-039 

U.14U_o.o30 

2 070 + 0.125 
z.y/y_o 190 

0.92/ 82 

(Mrk 335) 

2 

n 1 qo"I"0.055 

u.ioz_Q Q55 

U.UZ4_o 024 

- 

9 m 7-+0.078 
z.ui i _0.076 

U.lzu_o.oi3 

9 7Qr\+0.096 

Z. <^oU_0.084 

1 .01/120 


3 


0 Q74+0-053 

U.U^4_o 053 

- 

2 145+0.080 
z.iq:0_o.o77 

0 140'*"®'^^^ 
U.l4U_o.oi7 

2-9211;5+^ 

0.89/120 


4 

0.077l°;“g 

0 Q 7 n +0 061 

U.um _0 061 

- 

r\ i'7>7+0.093 

z.ii 1 _0.090 

0.l60l°;“® 

2.896l“;i®? 

1.05/112 

0006810101 

1 

U.iU/ _Q 

u.izi_o 046 

U.UUD_o 036 

9 4.79 + 0-037 

0_Q Q^Q 

0.189l°;°l2 

2 401 +0-085 

Z.4Ul_o.o76 

0.93/137 

(Akn 564) 

2 

0 qqq+ 0-048 
U.Uyo_Q 048 

n 1 o' 7~I"0.054 
U.iZ/_o 054 

0 041 +0-044 
U.U41_0.o44 

9 4QR+0 045 

Z.4yD_o.o44 

0 1Q7+0-017 

u.iy/'_o.oi5 

9 9 QQ+O.O 93 

z.zoy_o 083 

0.99/130 


3 

U.1^4_o 052 

0 1 qq+0059 

u.iyy_o 059 

U.U40_o 047 

9 401 +0.048 
Z.4tyi_o 047 

0 204"*'®'^^® 
U.ZU4_o 016 

c) i'7'7+0.092 
^• 1 ' '-0.082 

1.05/126 






